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OPTIMALITY STUDY OF A GUST ALLEVIATION SYSTEM FOR LIGHT 


WING-LOADING STOL AIRCRAFT 
Masaki Komoda* 

Ames Research Center 


SUMMARY 


An analytical study was made of an optimal gust alleviation system that 
employs a vertical gust sensor mounted forward of an aircraft *s center of 
gravity. Frequency domain optimization techniques were employed to synthesize 
the optimal filters that process the corrective signals to the flaps and ele- 
vator actuators. Special attention was given to evaluating the effectiveness 
of lead time, that is, the time by which relative wind sensor information 
should lead the actual encounter of the gust. The resulting filter is 
expressed as an implicit function of the prescribed control cost. A numerical 
example for a light wing- loading STOL aircraft is included in which the opti- 
mal trade-off between performance and control cost is systematically studied. 


INTRODUCTION 


The need for gust alleviation in reducing peak loads and improving ride 
quality has been apparent since the early days of aviation, and numerous 
studies of alleviation techniques have been conducted. A review of a number 
of these efforts and of the specific problems associated with gust alleviation 
is given in reference 1. In the last few years, the need to improve ride qual- 
ity for light wing-loading STOL aircraft in short-haul service has led to 
increased interest in gust alleviation systems. Analysis has shown that the 
severe disturbances in the longitudinal mode arising from vertical gusts can 
be effectively alleviated through the use of direct-lift devices, such as 
flaps, as the primary control (refs. 2-5), and that systems that employ rela- 
tive wind sensors in conjunction with inertial sensors (ref. 6) are more 
attractive than other types. The use of a relative wind sensor permits the 
gust velocity component to be extracted as a disturbance signal so that a feed- 
forward loop can be constructed that does not interfere with the pilot’s con- 
trol of the aircraft. This arrangement also permits the use of gust signals 
with positive lead time (see fig. 1). 

In most analyses of gust alleviation systems, the filter that lies 
between the sensors and the control actuators has been assumed to be of fixed 
shape and optima have been defined through variation of system parameters such 
as filter gains. Some alleviation system studies have been made in which the 
filter shape is not assumed, but is optimized using frequency domain analysis 
(refs. 7 and 8). In any practical system, control cost, such as surface 

*NRC Postdoctoral Research Associate at Ames Research Center. Presently, 
National Aerospace Laboratory, Mitaka, Tokyo, Japan 



deflection or deflection rate, has some limit. Since limited cost results in 
limited performance, it is of interest to examine how the use of an optimally 
filtered lead signal improves system performance. 

The present analysis extends the frequency domain optimization technique 
to define optima for systems with a feed-forward loop with lead time provided 
by a relative wind sensor located forward of the aircraft’s center of gravity. 
A numerical example is given, for a light wing-loading STOL aircraft, that 
illustrates the effectiveness of lead time in reducing control cost for sys- 
tems that use flaps alone, elevator alone, and combinations of flap and 
elevator. 


DESCRIPTION OF SYSTEM 


Major Assumptions 

The objective of this analysis is to improve aircraft ride quality. 
Vertical acceleration disturbances within a frequency band of approximately 
0.2 to 20 Hz are of particular importance from a ride-quality standpoint 
(ref. 9) . Analyses covering this frequency band require the inclusion of 
unsteady aerodynamic effects and elastic modes in the state equation. These 
effects are omitted in this analysis for reasons of simplicity. Even for large 
and flexible aircraft, however, the rigid body modes make a fundamental con- 
tribution to the total vertical acceleration disturbance; therefore this study 
should give a good indication of the possible ride-quality improvements. 

Spanwise variation of the vertical gusts is also ignored in this analysis. 
This may have a significant effect on the system performance (ref. 10) and 
should be considered in future extensions of this work. In addition, future 
studies should include the synthesis of the optimal system which gives the 
best estimate of the equivalent one-dimensional gust component under the 
presence of measurement noise. 

Airplane Dynamics and Turbulence Modeling 

Linear perturbation equations are assumed to describe the rigid body 
response of the airplane both as to control surface deflection and atmospheric 
turbulence. Considering a particular realization of turbulence which is 
temporarily assumed to vanish outside a long but finite time interval \t\ <_T/2 

x(a) = Fa(a)6(s) + 6^^(a)w^(a) (1) 

is obtained after bilateral Laplace transform using complex frequency a, where 
x = (w,a,0)^ is state, 5 = (6j?, 6^)^ is control surface deflection, = (Ug,Wgf 

is longitudinal mode disturbance, and Pa and are corresponding matrix 
transfer functions with consistent dimensions. Appendix A summarizes the 
expressions for F^^s) and 


2 


r 


The spacewise wave form of a frozen turbulence is approximated at each 
instant by a "linear representation." That is, the velocity component itself 
and only the linear part of spacewise variations of turbulence, both evaluated 
at CG, are considered. Thus, the rigid body aerodynamic derivatives are 
consistent in constructing . 

Since rigid body response is the primary concern in this study, no 
flexible mode is considered. By the same token, no unsteady aerodynamic 
effects are included. However, by incorporating the lift growth functions 
such as Wagner’s and Kiissner’s, into Pa (s) and (7^Cs), the refinement necessary 
to consider unsteady effects could be readily made. 

The disturbance power spectrum is given by 

'i’(w) = 7 J ^ fr/2 dx (2) 

where to is the angular frequency, is a realization of turbulence, and 

where E stands for the ensemble average of turbulence. Although the existence 
of cross power spectra and their effect upon airplane response are suggested, 
especially near the ground (ref. 11], only diagonal elements of ^>(co) are 
retained. A one-dimensional Dryden model of continuous turbulence is used in 
the analysis. The one-sided power spectra are given by 




7T 


diag 



2 

O S 


3(0^ + (0^^ 

(to^ + 


(3) 


where gust frequency, to^ = U/L^ and rms intensity, Og, are the primary param- 
eters. Excluding off-diagonal terms in $ (w] simplifies the optimization proc- 
ess because this separates the effects of Ug and w^y and the results thereof 

are additive. Since a similar formulation is valid for the Ug component, only 
Wq is considered in the following formulation. The term G(x{s)}rig{s') in equa- 
tion (1] is thus replaced by g^{s')Wg{s') y where g^(s] is the column of 6^^(s]. 

By ignoring unsteady lift growth effects and using a linear representa- 
tion for turbulence with a one-dimensional power spectrum, a trivial diffi- 
culty is introduced in the construction of Fa and if a conventional 

approximation e~^ = (1 - -s] is used for the down-wash and gust-wash lags (see 
eq. (A3)). These assumptions lead to an erroneously high response in the high 
frequency range. Actually, high frequency spectra of (s) and Gais^Vgis) 

are attenuated by both unsteady effects and spanwise variations of turbulence. 
To avoid this difficulty, while keeping a reasonable simplicity in the optimi- 
zation procedure, the first-order Fade expansion, e~^ = {-z + 2)/(s +2), is 
used to approximate these lags. 


3 


I 





Under the assumptions stated above, the transfer function matrices 
and Ga^s) have a common denominator (a) , that is the open loop character- 
istic equation. A stable open loop is assumed. Denoting by and G^y^^s) 

the numerator matrices of and respectively. 


^a(s) 




Gais) 


da^s^ 


(4) 


each entry of Fan and Gq^^ is a polynomial in s with an order not greater than 
that of dais'). 


Output Variable and Performance Index 

Suppose a scalar output variable y{s) is taken as the basis for assessing 
ride qualities, where y {s) is defined by 

z/(s) = h^(s)x(s) 

= f^(s)6(s) + g{s}Wg(^s) (5) 


with 


f^(s) 4 h^(s3^’^(s) 
i [/(s),e(s)] 

d(x (g) 


g{s) i h^(3)g^(s) 

Gn Cg) 

“ ^aCs) 


/ 


( 6 ) 


The output-defining vector h(a) is a polynomial in s. For example, if y {s) is 
the normal acceleration factor at a fuselage station (positive when 

rearward) , then 


h^(s) 



(7) 


In the case when y{s) is not a simple dynamic quantity but is, as an 
example, a bandpass filtered one, h(s) might be a rational polynomial in s. 
Additional poles would hence be included in f^(a) and gis). 

Since reasonable grounds for choosing otherwise ^ assessing ride 
qualities are lacking, the simple mean square value, is chosen as a per- 

formance index, Jy. Using the general transfer function y is) /Wg{s) y the 
performance index is then 
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• oo 

As is well known, the mean square, or equivalent, the rms, value of a random 
process can be related to other measures of the process such as the probability 
of exceeding a given level, even if the process is not Gaussian (ref. 12). In 
any case, the vector h(s) must be so chosen that Jy is well defined by 
equation (8) . 


2 

^ (a))do) (8) 


Control Surface Dynamics and Cost Index 
The control surface dynamics are generally written as 


6(s) = C(s)u(s) 


C9) 


T 

where u = (,Ujy,u^) is a hypothetically inertialess control input to the actu- 
ating system, and C(s) is the corresponding transfer function matrix. It 
would be reasonable enough to define the cost indices and by 


" 4/e 




K/e 


w 


9 


(jw) 




( 10 ) 


where and are the mean square input values of the flap and elevator 
actuating system and (s) (s) are corresponding transfer functions. 

Instead of trying to model actual systems precisely, a first-order 
decoupled dynamic system is assumed for C(s) in the analysis. Thus, 




diag 



s 



( 11 ) 


where and are non-negative real, and distinct. Using equations (11) and 
(9) , equation (10) becomes 


J. 


u 


//« 



JO) + Pjyg 

^f/e ^9 


(jw) 


2 




-2 


'f/e 


( 12 ) 


which indicates that the cost index defined above affords a good measure of 
the mean square values of control surface deflection and its time rate. 
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Control Law 


One of the basic assumptions of this study is that an instantaneous 
measurement of disturbance is available with some lead time The 

simplest way to measure the ^disturbance before it actually hits the lifting 
surfaces would be a combination of a relative wind sensor, mounted ahead of 
the wing, and inertial sensors that generate signals due to non-gust relative 
wind components. Although it is necessary to examine carefully the possible 
adverse effect due to measurement noise, more sophisticated versions, such as 
a predicting filter, could be utilized to increase the lead time 

Let it be assumed that there exists an optimal filter k(s) such that a 
control law 


u(s) = (13) 

minimizes the performance index Jy subject to the prescribed cost indices 
and Only asymptotic stability and causality are required for k(a) at 

this stage. Successive substitutions of equation (13) into (9), and of equa- 
tion (9) into (5), give the system transfer function with feed-forward loop 
closed. 


y(s) = g^is^Wgis) 

= [g{s) - f^(s)C(s)k(s)e'^9’®]w^(s) (14) 

Equations (13) and (14) are valid for each realization . 


OPTIMIZATION 


Augmented Performance Index 


Following standard procedures, an augmented performance index J is 
defined by 

J - Jy ^ X/Juf + ^e^Jue 05 ) 


where and X^ are real multipliers. After the substitutions of equa- 
tions (14) and (13) into equations (8) and (10), respectively. 




x/j 

J uf & ue 


1 r 

1 r“ 'p 

J Tr j Ak(-ja))k^(jw)$w (o))<f(o 

^ — oo y 


(16a) 


(16b) 
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where 


A = diagiX^^ (17) 

Using equations (16a) and (16b), equation (15) is expressed as a functional of 
the unknown filter k(jo)). 


Weiner-Hopf Equation 

Instead of solving for k(a), it is more convenient to define and solve 
for a vector 


v(s) = kis)Wg^[s') (18) 

where (s') is a factored form of ^^j_(a)) in the complex s-plane, that is. 





(19) 


where, from equation (3) , 


W^g(s) k Og^ 


S + 0J g//3 
(s + w^)^ 



(19a) 


which is usually referred to as the equivalent deterministic disturbance. As 
is known, the original (stochastic) problem of finding k(a) is equivalent to the 
deterministic problem of finding v(s) with respect to the same J. Once v(a) 
is found, k(s) is readily obtained by 

k(s) = v(s)w^^(s) (20) 


As derived in appendix B, the Wiener-Hopf (W-H) equation to solve is 


[z(s)]^ £ [/l(s)v(s) - b(s)e~’^^®]_^ = 0 (21) 

where 

4(s) 4 A + C'(-8)f(-s)f^(s)C(s) (22a) 

b(s) 4 Ci-3)U-s)g{s)Wg^{s) (22b) 


and where [ indicates that part of [ ] which is analytic in the RH s-plane. 
Since the matrix i4(s) is real and positive on the imaginary axis of the 
a-plane, solution v(s) gives a minimum of J. 
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Solution v(s) 

As shown in appendix B, the equivalent optimal control v(a) is given by 
(ref. 13) 


v(s) = p-i(s)[[p2’(-s)]-lb(s)e'V]_^ (23) 

where P(s) is a factored matrix of A(s) ^ such that P(s), as well as P~^(a), is 
analytical in the RH a-plane. Contrary to the apparent simplicity of equa- 
tion (23) , it is not easy to find the factored form P(a) except when only a 
single controller is assumed. Here, ^4(s), and hence P(a) and v(a), degenerate 
into a scalar rational polynomial. The case of a single controller is 
summarized in appendix B. 

An alternate method of finding v(s), including the case with Xg > 0, is 
summarized in appendix C. As derived there, the equivalent optimal control 
v(a) can be decomposed into two parts, 

v(s) = v^(s)e”’^9'®+ v^(s;t^) (24) 

where v^(s) is the basic part of v(s) when = 0, and Vp(s;xg) is the 
additional part when xg > 0, 

To study the properties of v(a), let us consider the poles of the system 
with feed- forward loop closed. The fundamental characteristics of the system 
are governed by a set of system poles s - k = 1, . . ., ky^, which are LH 

zeroes of det[A{s)'\ ^ that is, with Re{syid) > 0, 

det[Ai-Suk)] = 0; k = I, (25) 

With a scalar output variable and with the first-order control surface 
dynamics, the number of system poles is given by 

ky - ia (number of controllers) (26) 

where is the order of the least common denominator of f(a). Since -4(a) is 
Hermitian, equation (25) indicates that s - Syj^\ k = 1, . . ., ky, is also a 
zero of det[i4(s)]. Let us define a polynomial 

duis) Ml (s + s^-0 (27) 

fe=l 


for later use. Poles Sy]^ are neither dependent on the nature of the disturb- 
ance nor on Xgy but are dependent on the prescribed values of cost via 
multipliers Pf/^ and Some pertinent properties of v^(a) and Vp{s\Xg) 

will follow. 
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Basic Part v^(a) 


From equation (C28a) , 


v^(s) =^7-l(s) 


(s + Oig)^ • duis) 


(28) 


where the numerator polynomial p^^(s) is at least two orders less than that of 
(s + , so that v^(s) is strictly proper. This part represents the 

equivalent optimal control when the gust sensor is located at the CG, or equiv- 
alently Tg = 0, and is invariant whatever 0) is. The corresponding part 

k^(a) of optimal filter k(a) is given by 


k^Cs) = C-l(s) 


Pon(®) -I r ^ 

(s + 0)^)2d^(s) 


(29) 


= C. 


Cy 


s + 


/3 


E 

k=i 




^ ^uk 


(29a) 


Equation (29a) is the partial fraction form of equation (29) after extracting 
its constant part c^. Equation (29a) indicates that k^(s) is composed of a 
train of first-order (including Cq term) and second-order (for complex conju- 
gate low-pass filters. Each gain Cq, c^, and has been so determined 

that the best compromise between performance and costs results. Figure 2a 
shows schematically the structure of k^(a)e''^9'^, in the corresponding impul- 
sive response function ^ [k^ (s)e“'^9'^] representation. Each element of filter 
Y^{s')e~'^9^ has an impulsive response that vanishes when '^g ^c>r > 0. 

The corresponding system transfer function ~ with the 

feed- forward loop closed is given by (see eq. (C33)) 


9uo (s) 




^uon 


where the order of is equal to that of (s + ly^g/ZZ') 

\xjy • 0 . 


dy{s') when 


(50) 


Additional Part Vp(s;T^) 

The part Vp(s;T^) is that part of v(s) which makes maximum utilization of 
the information about the turbulence that has passed the wind sensor station 
but has not yet reached the CG. From equation (C28b) , this part is best 
expressed by 
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(31) 


hj V 


VpCSlTg) =C-1(S)^ 


- e~^9^uk ^ Yp^(O)^'''^^ - Yp;^(V 


^ok 


-s + Suk 


® + ^uk 


where 'Ypk^.'^g^ coefficients dependent on Tg and where Ypfc(O) = Yp;;/(T^ =0). 

The coefficients are invariant with t^. Obviously, Vp(s;T^), and hence 
the corresponding part kp(s;Tp), are not a finite order system. This is so 
because of the lag terms that are present in equation (31) . Some con- 

sideration of a rational polynomial approximation is contained in a later sec- 
tion. Also, clearly Vp(s;xp) tends to vanish when -> 0, thus leaving only 

the Vq(s) portion of equation (24). It must be noted that Vp(s;xp) is finite 
and analytic in the RH s-plane, even at s = k = 1, . . ., ku- The 

counterpart kp(s;xp) of k(s) is written as 


kp (s ; Tg) 


C'-l(s)[C(s)Vp(s;Xp)]w'^(s) 

P- -x^s 7 r '1 1 ^ - e 

[lk(0-)e 9 -^^(Tp)]+r^ 

+ ^ ~'^k^'^9^ ^ s^(0)e ^ - s^(Xp) 

s + aj^//3 ^ ^uk 



C32) 


(32a) 


where gain vectors ^^(r^), readily obtained as 

linear combinations of 3^^ and equation (31). (See eq. (C37) for 

details.) Thus, kp(s;r^) as seen in (32a) constitutes another train of linear 
filters. Figure 2b depicts the impulse response of each term of equation (32a). 

As seen, the impulse responses vanish for x < 0, and hence k^(s;x ) is again 

^ y 

causal. For x > Tg, kp(s;Xg,) is exponentially stable. 


Optimal Performance and Costs 

As a result of the above, one has the optimal filter of the form 

k(s) = k^(s)e"V + kp(s;Xp) (33) 

When equation (33) is substituted into equations (16a) and (16b) using equa- 
tion (14), the performance Jy and costs 9ufle> which are optimal in the stated 
sense, are obtained. Integrals in equations (16a) and (16b) are evaluated by 
summing up all residues in the LH or RH s-plane. The resulting Jy and Jy^fjQ 
are implicitly dependent upon the assumed value of the multipliers and 

the actuator time constants yy/g (see eq. (11)). One can then see the best 
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trade-off between available performance and required costs with the ]i and X 
as parameters. 


Suboptimal Filter k(s) 

In physical implementation of equation (33), some difficulties may arise, 
especially in the r^-terms in equation (32a). To cope with these difficulties, 
a rational polynomial approximation of k(s) seems preferable. A simple way to 
do this is to approximate the e ^9'^ terms, as well as the terms by a 

Fade expansion (ref. 14) . It is expected that the higher the order of the 
polynomial approximation, the better the approximation will be (ref. 15) . 

Just for simplicity, the first-order expansion formula 




-z 


-z + a 
z + a " 


z : complex 


(34) 


is assumed hereinafter, with a real and positive constant a. Substituting 

= (-S + a/T^)/(s + a/x^) etc. into equation (33) and executing the nec- 
essary manipulations, one obtains a suboptimal filter of the form 


-s + 


k(s) = k^(a) 


s + ■ 


a 


^ + icp(s;T^) 


'^9 


where k^(s) is given by equation (29a) and where 




r(T^) t(0 


§7 (t ) 

k'' g-’ 

• ® + ^uk 


"9 


9-' 


s + 


a 

'^9 


/3 


(35) 


(36) 


The gain vectors ^(x^), r(x^), t(x^) and s^(x^) are uniquely determined as a 
linear combination of and yp'ki'^g) > where Yp?<(x^) stands for the Tp^('fp) 
which is solved with the Fade expansion of equation (34) in equations (C25b) 
and (C25c) . It is not difficult to see that the suboptimal filter k(a) is the 
optimal filter of the W-H equation (eq. (21)) to which the same Fade expansion 
has been incorporated at the beginning. Corresponding suboptimal performance 
Jy and suboptimal costs are obtained by equations (16a) and (16b) as 

before. 


NUMERICAL EXAMFLE AND DISCUSSION 


A typical light wing-loading STOL plane is chosen in this example. The 
pertinent numerical data are summarized in table 1. Since we are studying the 
nature of the optimal filter and the resulting gust alleviation, both with and 
without lead time x^, only one flight condition (£/ = 41.45 m/sec) and one 
turbulence scale length (L = 305 m) are assumed. The necessary aerodynamic 
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derivatives are summarized in table 2. The normal acceleration factor at 
three fuselage stations (1^ = -2.5 m, 0.0 m, and 2.5 m correspond approxi- 
mately to the cockpit, CG, and aft cabin stations, respectively) is taken as 
the output variable j/(a) (see fig. 1(a)). However, primary interest is placed 
on the aft cabin because it is at this station that alleviation is usually the 
most difficult to achieve. 


TABLE 1.- AIRPLANE DATA AND FLIGHT CONDITION 


Mass, kg 4,987 

Wing area, m^ 39.02 

Mean aerodynamic chord, m 1.981 

Radius of gyration about pitch axis, m 2.576 

Tail moment arm, m 7.81 

True airspeed, m/sec 41.15 

Altitude, m 0 

Turbulence scale, m 304.8 

Flap setting, deg 20 


TABLE 2.- AERODYNAMIC DERIVATIVES^ 

^6 f 

— = 0.465 
ea 


Derivatives 

Cijf 

Cv 


Cm 

a( )/9u 
9( )/9a 
9( )/36f 
9C )/96g 

0.12 

-0.012 

0.12 

0.521 

0.344 

0.1214 

6.42 

2.58 

0.453 

-1.26 

-0.195 

-1.78 


— 

— 

— 

-7.97 

9( )/»(§) 

— 

— 

— 

-23.95 


^Definitions of nondimens ional derivatives 
in table 2 and of stability derivatives in 
appendix A are found in reference 18, where a 
conversion w = Ua is used. 
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Flap Control With Xg = 0 

Figure 3 shows optimal trade-off characteristics between and 

when Tg ~ 0 with the flap actuating system time constant 1/p^ a 

parameter. As is known, a flap system (or any equivalent direct- lift con- 
troller) is very effective in alleviating the response to vertical gusts when 
a reasonably quick response (or large and enough power (or large v/u^) are 

available in the flap actuating system. When the available rms input power 
^/^ increases (i:e., decreases), the corresponding system poles (^ = ~s^^) 

move as shown in figure 4. Ignoring one real pole, which corresponds to an 
artificially introduced downwash lag pole, two pairs of complex conjugate 
poles modify the basic short period and phugoid poles. Another real pole 
modifies the surface actuating system pole (a = -yj?) . 

In this particular example, and in most cases with a well designed direct- 
lift controller, the 6^-into-ng transfer function /(a) has only LH zeroes. 

J 2 9 

Hence, considering a limiting case Xj? 0 (with X^ “) , the system poles 
approach the open-loop zeroes (see eq. (Bll)), and an optimal control law 


<S/(s) _ ffnCg) 


WgCs) 


/„(s) 


(37) 


is obtained (see eq. (B13)). Equation (37) defines a stable filter w hich 
could perfectly cancel the disturbance (see eq. (B14)) with a finite when 




or with an infinite when < 


Practically, finite values of iij? and must be considered. As seen in 

figure 3, a larger time constant requires more rms input power to 

achieve a given performance. The same may be seen from a different point of 
view in figure 5 ( curv e for = 0) where the trade-off characteristics 


between and 


( curv e for 

, which are required to achieve = 0.03 g at the 


^ y Wli-UV-'ll X -L J. VX CX 1 i X V-/ V r ' 

cabin station when a = 2.1 m/ sec, are shown. Figure 6 shows the system trans- 
fer function 9'^^(a) with the feed- forw ard loop closed, corresponding to 
several combinations of and ^jy^ in figure 5. From figure 6, one can see 

the frequency-wise structure of the achieved alleviation. As seen, a smaller 
time constant 1/y^ provides more alleviation at the high frequencies, but less 

alleviation at the low frequencies. Figures 5 and 6 clearly show that a 
clever choice must be made from the ride-qualities standpoint to achieve a 
proper comp romi se be twee n the frequency band to be alleviated and the 


available ^jy^ and 




Flap Control With > 0 


As indicated in app endix B (see eq. (B13)), lead time x^(>0) is 
unnecessary when enough is available. In other words, if the flap 
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surface can be actuated quickly enough and the effect thereof is of minimum 
phase lag, then the corr ectiv e control need not be initiated before the dis- 
turbance arrives. When is finite however, lead time Xg reduces the 

required cost to some extent. Examples are shown in figure 7. Although it is 
unrealistic, the limiting case -> » is included to indicate how lead time Xg 
affects gust alleviation. When l/v^f is zero (fig. 7(a)), only a negligible 


cost saving is available with a realistic value of = 0.1 sec. When 1/Vj? is 

0.5 sec (fig. 7(b)) some improvement is realized by small lead values 

{xg =0.1 sec). These cost savings are also compared in figure 5 with l/Mf 

and Xg as parameters. 


Generally, the merit of lead time Xg is only realized when the actuating 
system response is slow. This fact would be further emphasized when a sub- 
optimal filter is substituted for the optimal one. The broken lines in fig- 
ure 7 show the performance corresponding to a rational polynomial filter as 
given by equation (35) with a = 2. In these examples, apparently the disad- 
vantage introduced by rational polynomial approximation destroys the merit of 
positive Xg. This shows that a higher-order Fade approximation would be 
necessary. 


Elevator Control With x^ = 0 

Figure 8 shows the optimal trade-off between and /wg^/a^with 

x^ = 0 at the three fuselage stations. As is well known, the elevator is not 
an effective gust alleviator, especially at the aft cabin station. This is 
due to non-minimum phase lag characteristics of the Sg-into-n^ transfer func- 
tion e(s). Figure 9 shows the system poles s = with /Wg2 (or equi va- 

lently with Xg^) as a parameter at the cabin station. When allowable 
increases, similar to the flap case, system poles start from open loop poles 
and approach the zeroes of e{s') . When there is a non-minimum phase zero, such 
as shown in figure 9, one of the system poles corresponding to the short 
period mode terminates at the image in the LH s-plane of the RH zero. This is 
what is indicated in equation (B16) . The limiting case 0 (with °°) 

gives an optimal control law (see eq. (B18)) 


Wgis) 


1_ _ 9n* (s) 



(38) 


where (s) is with replacement of its RH zero factor (-a + by its 

image factor (s + , and where p^*(s) is a polynomial defined by equa- 

tion (B20) . Equation (38) defines another stable filter, but, in this case, 
perfect alleviation is not obtained. In fact, the corresponding system 
response is given by (see eq. (B21)) 
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( 39 ) 



In figure 10, which shows system response with as a parameter, one 

can see the frequency-wise structure of the alle viation. It is clear from the 
figure that the desired alleviation in rms is obtained only at the sacri- 

fice of amplified response in the low- and high-frequency ranges. Finally, 
the best available performance corresponding to the limiting case 0 is 

obtained from equations (16a) and (39) , 



which is achieved with infinite cost /wg2 , Using equations (40) and (19a), 
equation (41) is rewritten in a more practical form: 




(Sg + 0)gr) ^ 




(41a) 


The inefficiency of elevator control as a gust alleviator depends on the 
fuselage station because the location of the RH zero s = Sg is dependent upon 
the station as well as on the flight conditions. The limiting performance of 
equation (41) without lead time in gust sensing is plotted in the ~ ^ 

curve in figure 11 against fuselage station. At the fuselage station aft of 
the center of percussion positive rearward from CG) , the 

ride quality becomes worse. At the stations before Icc^ there exist no RH 
zeroes in e^(s), hen ce t heoretically perfect alleviation is achievable with 
enough input power /ug2 as in the case of flap control. Using the short 
period mode approximation, the marginal point (center of percussion) is 

given by 



Elevator Control With > 0 

Let us consider the same limiting performance when Ag^ ^ 0 with > 0 
at fuselage stations la > lac* derived in appendix B, the additional part 

of the filter adds the corrective control 
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ey* (s) 


( 43 ) 


<Se(s) 


e"V - _.-l 


•s + s. 




and this results in a total system response (see eq. (B23)) o£ 




(44) 


and the limiting performance as /w^2 -> oo of 


^ 9 ^ 
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(45) 


Although equations (44) and (45) are only realized with an infinite control 
input they suggest the possible improvement that lead time gives, an 

improvement that can be significant especially when Zq is large. The essen- 
sial role of 6^(s)/zJ^(s) in equation (43) is that the RH zero in e(s) is can- 
celled out by the same factor in the filter denominator when Oyipis) is 
considered (see eq. (B22)). However, when a finite cost (with > 0) 

is prescribed, such a pole-zero cancellation occurs only partially. The 
effect of the additional part of the filter is shown in figure 11 with > 0 
as a parameter. Figures 12(a) through (c) show the optimal trade-off between 
and when > 0 at the three fuselage stations. The ride 

improvement is considerable when ^ 0.1 sec at the aft cabin. As is pre- 
dicted by equation (44), the alleviation due to a filter which includes the 
additional lead time part covers a wide frequency band. This is shown in fig- 
ure 13 where system response 0'^(a) with and without lead time is compared. 

The performance with suboptimal filters as given by equation (35) shows 
no significant performance degradation with a = 2. This suboptimal filter is 
compared with the optimal filter in figure 13 for the case of A^ = 0.1 and 
= 0.1 sec. 


Two Controllers 

It has been shown that the flap control is effective for gust alleviation, 
if enough control power is provided, and that the effectiveness of elevator 
control as a gust suppressor is improved by feeding the disturbance signal in 
with lead time. Now let us consider how much the required flap system cost 
could be reduced by incorporating the elevator as an auxiliary controll er. 
Figures 14(a) and 14(b) show optimal trade-off characteristics between 

/u^^logy and with Xg as a parameter. These figures are results of the 

exact filter of equation (33), but the results using the suboptimal filter of 
equation (35) with a = 2 differ only indiscernibly from those shown. This is 


16 




constant 1/^^. (See suboptimal performance in figs. 7(a) and 7(b).) Similar 

to the single controller case, the elevator shows some effectiveness at the 
cockpit station even when Tg - 0, but it makes little contribution at the aft 
cabin station. When '^g the situation is improved considerably at the 

cabin station. 


Hereinafter, attention will be focused on the aft cabin station, where 
non-minimum phase lag characteristics hamper the effectiveness of elevator 
control. The cost savings achieved in the flap system by incorpo rati ng the 
elevator, is shown in figure 15. In the figure, required costs and 

JU to obtain = 0.03 g for = 2.1 m/sec are plotted with rms elevator 

power and flap time constant l/py* as parameters for '^g - ^ sec and 


0.1 sec. These results are based upon the use of a suboptimal filter as in 
equation (35) , and hence the adverse effects of rational polynomial approxima- 
tion appear when 1/Py^ is small with Tg > 0. Except for this fact, it is seen 

that much greater cost savings can be achieved in the flap system by elevator 
control when > 0. As to the frequency-wise structure of the alleviated 

system, several cases are compared in figure 16 in their system transfer func- 
tion g^{s') , The figure indicates that the resulting alleviation by two con- 
trollers is the blended one of the single controller cases of figures 6 and 
13. This is shown by the fact that the corresponding system poles s = 
are located between the optimal pole locations in the single controller cases. 
This is seen in figure 17 where the system poles corresponding to the short 
period mode are shown. Locations of the other poles are almost the same as 
those in the single controller cases. 


Limiting Cases 

As h as b een seen in the above, if one is concerned with a given criteria 
such as ^^lOg ~ 0.03 g/2.1 m/sec at a particular fuselage station, there is 

a family of required co sts (yj:*, /li^) and (y^ , , or equivalently 

(/6j?^,/l^) and as illustrated in figure 18. To obtain a rough 

idea of system performance, including the power spectrum characteristics and 
required costs, two limiting cases are considered. One is the case of flaps 
only (Xg^ ^ "^) ; it corresponds to point F in figure 18 and has been described 
before. The other is the case assuming an infinite elevator control. It is 
indicated by point E in figure 18 and is characterized by 0 while 

retaining XjP- finite. A process similar to the one used in appendix B 
(limiting case A^ ^ 0) leads to somewhat simple formulas. Inspection of 
equation (C16) gives the limiting characteristics 

at ^ c s + ]Sf y^ 

vduCs) Xf °° • ^ %*(«) ; ~ 06) 

J ‘-00 '^e 
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where eyip is the leading term coefficient of and where has the 

same definition that is used in equation (38). Also, equation (C19) yields 




s + Pg 


g-nis) 


\_ Ve ' -s + ■ s + ■ e„(s)“ ' 


( 47 ) 


Denoting by (s) and y®(s), and p-^(s) and p®(s) the flap and elevator compo- 
nents of v(s) and p(s), respectively, equation (47) indicates (see eq. (C9)) 


/ (s) 






P-' (S) 


u®(s) = 


s + y. 


t Ooo 

IFT) 


s + <?„ 




w^g(s)e'V -pe(s)j 


(48a) 


(48b) 


In equation (48a), p-^ (s) has no poles other than du{s) zeroes, which comprise 
the LH zeroes of e^*(s) as well as s = -py- and s = -g„. Evaluating equa- 
tion (C25a) at these poles, it can be shown that all the unknown coefficients 
in p-^(s) vanish except for the ones, ^^d , which are the gains 

corresponding to the poles s = -yy> and s = -Sg, respectively. Further, equa- 


tion (C25c) requires Yp^(.'^g) = 0, so that one can write 


)/(s) . 




s + y 


f 


s + 


(S + Py.) (S + 3g) 


(49) 


Now considering equation (47) again, the only possible RH pole in p®(s) is 
s = 3_, which comes out of e^(s) . This implies that the RH members of equa- 
tions (C25b) and (C25c) are proportional to e S' and hence that all the 
unknown coefficients (t^) Yp^'^g') given by Yp(0) • e ~^g^^ and 

Yp(0)e”^S^®, respectively. Thus, substituting equation (49) into equa- 
tion (48a) , one has 

f 

yj(s) = - Ve (50) 




s + z. 


which gives the required cost in the flap system 
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(51) 




(0) • 




The elevator component v^(s) has a similar structure to that of equation (B18), 
except that it includes an additional pole s = -yj?? and that both and 

are now dependent upon • The corresponding*^ system transfer function of 

equation (14) is written as 

Fv^ (afj 

gM = ^(s) - [/(s)e(s)]C(s) \w~Us^e^g^ (52) 

[y®(s)J 


Substitution of equations (48a) and (48b) into equation (52) gives a limit 
when a ^ 

cx» 


g-j^Cs) -> [/‘(s)p-^(s) + e(s)p®(s)J w^^(s)e’^9'^ 


(53) 


Taking the requirements of equations (C25a) through (C25c) into account, one 
can expand the RH [ ] of equation (53) into partial fractions and find that 




(54) 


where the system gain pg^ is given by 


Xj- zj- - y 2 


ef 


. IL 


— y{ fo) 


fi^e) P« 


(55) 


Equation (54), a very similar expression to equation (44), which is the limit- 
ing case of elevator control, indicates that the output y{s) = has a 

power spectrum 




ef 


(o2 + 


e ^'^g^e 


(56) 


which is to be compared with the output power spectrum of open loop 


[3>ng(“)]o = 9'(-iw)^u^(“)S'(» 
Finally, equation (56) gives a limiting performance 


(57) 
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Since the filter gain Yp^CO), and hence the system gain is implicitly 

dependent on the assumed value of Xj? as well as on equations (51) and (58) 

do not afford a pair of closed form solutions of the flap cost and performance. 
However, it may be noted that the algorithm required to solve for the unknown 
gains is greatly simplified in this limiting case. 

An example of the reduction in the flap cost available by incorporating 
elevator control to achieve a given performance has been shown in figure 15. 

It is seen in figure 15 that a much greater cost reduction is realized in the 
flap system when the elevator control is used together with a positive lead 
time (e.g., Tg = 0.1 sec) than when it is used without lead ti me ( t ^ = 0 sec). 
Figure 1 9 co mpares the output power spectra corresponding to =0 (point F) 

and to oo (point E) with that of open loop. When the loop is open, two 

peaks appear in corresponding to the phugoid and the short period 

modes. The turbulence spectrum ^.nd the transfer function g{s^ cut off 

the output spectrum at the high- and low-frequency ranges, respectively. The 
same feature exists in the flap case (point F) where optimally relocated poles 
and zeroes result in a gust alleviation. In the two controller case (point E), 
the infinite elevator cost eliminates all the inherent poles and zeroes except 
one image pole s = -Sg. Thus the output spectrum (to) is not attenuated in 
its low frequency range. 


CONCLUSIONS 


An analytical study was made of an optimal gust alleviation system that 
uses vertical gust sensors mounted ahead of the CG. Frequency domain optimiza- 
tion techniques were used to synthesize the optimal filters that produce the 
corrective control signals into the flap and elevator actuators. A special 
emphasis was placed upon the effectiveness of the time by which sensor informa- 
tion leads the actual encounter of the gust in reducing the rms value of verti- 
cal acceleration caused by the gust. The resulting filter, a lagging filter 
when the lead time is positive, is expressed as an implicit function of the 
prescribed control costs; that is, prescribed rms values of surface deflection 
and its time rate. Using this method, the optimal trade-off between system 
performance and required costs was systematically studied. A numerical 
example, which considers a light wing-loading airplane in the approach condi- 
tion, is included. The conclusions derived from the study are that optimal 
feed-forward of gust signal with positive lead time results in the following. 

1. It improves the cost-effectiveness of flap control, especially when 
the flap actuating system has a large time constant. 

2. It improves cost-effectiveness of the elevator control significantly, 
especially at the aft cabin stations where otherwise non-minimum phase lag 



characteristics hamper the elevator’s effectiveness. Elevator control both 
with and without flap control decreases the rms vertical acceleration by a 
factor times that with the zero lead time = 0 case, if enough 

elevator control power is available. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, August 19, 1975 
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APPENDIX A 


STATE EQUATIONS 


The longitudinal equations of motion defining the response of an aircraft 
to control inputs and disturbance are written as 



-^a 

g 


-Zu 

^ - ^a 

-s 

a(s) 

.-Mu 

-iM^s^M^) 

sis - Mgl 

J Cs)_ 
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\J^u 
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Ug{s') 

Wg{s') 


where the linear representation of the vertical gust 


a - ^ 
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• _ i ^ 
V u u 


_ 1. ’^^g 


9w, 


3a; 


U U 


CAI) 


(A2) 


is assumed. Also assumed in equation (AI) is the first-order expansion in 
downwash lag terms such as 


a(t) 




Ziji 


u 


ait) 


(A3a) 


or equivalently in the frequency domain 

j^l - e" = -^sa(s) 


(A3b) 
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Inspection of the moment equation in equation (Al) indicates that response 9 
due to 6 j? or Wg, or both, becomes very high in the high frequency range. This 
is of no*^ consequence in the conventional analysis of rigid body modes; however, 
if one is concerned with the rms value of a response evaluated over the fre- 
quency range [0,»), this unrealistic feature must be avoided. For this reason, 
the first-order Fade expansion + 2)/(s + 2) is used in this analysis 

so that equation (A3b) becomes 

W 

L s + 


Consequently, terms M^s , s, and - Mg)s in equation (Al) are multiplied 
by the factor (2U/l^')/[s + (211/1^)], 

By multiplying the inverse of the LH coefficient matrix in equation (Al) , 
explicit forms of F^(s') and are obtained. 
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APPENDIX B 


THE WIENER-HOPF EQUATION AND ITS SOLUTION 


Wiener-Hopf Equation 

Given is an augmented performance index J (eqs. (15), (I 6 a), and (16b), 


^ '^e^^ue 


; oo 

{ [g'(jo))z<;^^(ja)) - f^(ja))C(jcj) v(jcj)e'^“’^S'] 


x[w^g(-jw) 9 '(-jw) - v^(-jtj)C(-ja))f (-io 3 )e 


+v^(-jaj)Av(jo)) 


(Bl) 


Required is a stable and causal v(s) such that minimizes J. The first 

variation 6J due to an infinitely small variation 6 v(jcj) is given by 


= J [ 6 v^(-< 70 )) z (joj) + z^(-jo)) 6 v(ju) ] (B2) 

with 

z(joj) = [A + C(-J 0 ))f (-jaj)f^(jw)C(jw) ]v(joj) 

-C(-jco)f (-ja)) 0 '(ju)w^^(jto)e"'^"’^S'“ (B3) 


In order that 6c7 = 0 for any variation 6 v(s), which is analytic in the RH 
s-plane, z(s) must be analytic in the LH s-plane, see, for example, refer- 
ence 9; this is stated by equation (21). It must be noted that, if v(s) is 
analytic in the RH s-plane, then so is k(s) of equation (20) and vice versa. 


Solution v(s) 

Since the matrix A(s) in equation (22a) is real on the jw axis, conjugate, 
and Hermitian, and has a maximum rank of 2 (when Xf^ . 5 ^ 0) or 1 

(when Xjy^ or = 0), there exists a (2x2) or (1x1) matrix P(s) such that 
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I 


yl(s) = P^(-s)P(s) (B4) 

and that P(s) as well as its inverse P”^(s), is analytic in the RH s-plane 
(ref. 13) . Hence, the solution v(s) is given by the closed form 

v(s) = P“^ (s) [[P^(-s) (B5) 


Single Controller Case 

Let either /(s) or ^(s) of equation (6), and X 

and y be corresponding multipliers. Matrix A(^s) and vector b(s) of equa- 
tions (22a) and (22b) become scalars, and the factored form P(s) is readily 
given by 


vd^^Cs) 

~ dais) dais) 

where, using a shorthand notation d^(a) = s + y and a real constant v, 
'o‘^dui-s)duis) = \‘^dai-s)dai-s)dais)dais) + \i^fmi-s)f-^y^is) 


CB6) 


(B7) 


and where d^{s) has only LH zeroes s = ^ 1. 

After partial fraction expansion, the rational polynomial [P^(-s) ] ” ^b (s) is 
expanded as follows 


[P^C-s)r bCa)e 9 ' v ' ^ 


where the 
that 




•a' 








\ 

(3,1 

+ k^i ^ ^ukf 

(B8) 


^ ... term has poles of dais) and of Wg^is) . Knowing 


r _ f” / J_ p* 

L-® ^ ^ukL Jo ^ J-j 


“ ts 


-s + s 


uk 


e ds\dt 


)• 


''^g^ - e '^g^uk 


-s + s 


uk 


(B9) 
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the solution of equation (B5) is written as 


y(s) 


j If ^ 

Vv; y j L d^i-s)d^(.sy^9e J + 


-T^S 


ku 

k=i 


e~^9^ - e'^d^ uk J 
* ^uk j 


(BIO) 


After another partial fraction expansion, equation (BIO) is rearranged into 
the standard form in equations (24) , (28) and (31) . 


Limiting Case 0 


A special case where an infinite cost is allowed, is specified by letting 
X^ -> 0. If finis) has no RH zeroes, equation (B7) implies that the system 

poles approach either LH zeroes of /i„(s) or s = -o^ when X^ -> 0, so that 


vdj,(s) y/i^(s) 


s + e„ 




= 


ino 


(BID 


where fino is the leading term coefficient of /ln(s) . Equation (B8) 
simplifies to 


[P^(-s)]-lb(s)e"^5'® 


-s + c„ 


9-nis) ^ ^ 


TrrS 


(B12) 


Applying [ ]^. operator to equation (B12) while retaining X^ as finite and 
finally letting X^ 0, one can see that the contribution of c^/(-s + c„) , 
which is the only possibility yielding the additional part of the filter, 
vanishes even when > 0. Hence, the limiting case without the non-minimum 
phase zero is obtained as follows: 


u(e) 


y 


g + y 

S + (3oo 


9ni^^ 


(s)e~^9^ 


s + y 


Pn(s) , 
f„(s) ^9e 


(s)e~^9^ 


(B13) 

(B13a) 


which gives, in turn 


9^is) ^ 0 (B14) 

where ^^^(s) is the system response that is defined by equation (14). 


26 



Next, suppose has a RH real zero, s = z^, with Sg > 0. Putting 


/inCs) = (-S + 3 e)/i„(s) 


(B15) 


/j^(s) has LH zeroes only. This time, equation (B7) has a limit 


vcZy(s) ^ y(s + 


s + 


(B16) 


and hence. 


[P^(-s)]-ib(s)e '^9^ 


s + Ze 9n(.^) 


WaMe^'G^ (B17) 


-s -s + Sg da(.s) ■ 9e 


As a result, from equation (B5) , a limiting case with the non-minimum phase 
lag is obtained as follows: 


v(s) 


s + y 


9n ^ (^) 


L(S + + 03^)^ 




+ P 


da(s) 


^ [s + 


+ sg J ^ 


where 


Pe = 2sg[3-(s)w^g(s)] 


3=3.: 


|s| < « (B18) 


(B19) 


and where g'^*(s) is a polynomial of order such that 

(-3 + 3g)9'^*(s) = Og^is + oig//5)(s + 3g)9'^(3) - pg dais){_s + (B20) 

Equation (B18) , in contrast with equation (B13a) , comprises both the basic and 
the additional parts of the filter. The basic part in equation (B18) gives a 
system response 


&(s) - /(s) 


9n Cs) 


= — - — w^hs) (B21) 

(3 + 2 g)/*^( 3 ) (3 + 0)^)2 9e s + Sg 9e 


and the additional part in equation (B18) gives an additional response 
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d(2 (s) 

^ -/Cs)Pg (s + zs)f*y,is) ' 


(B22) 


-TgS 


e'^g^e -1 




-S + Zr 


Wg^iB^e g 


= - Cl - eV)w-'(s) ^^22a) 

s 

It must be noted that, in equation (B22) , the factor (-s + 3g) in the numera- 
tor of f(s) is cancelled out by the same factor in the filter denominator. 

This is the essential role of the additional part of the filter. From equa- 
tions (B21) and (B22a) , 

e^g%g^^(s) (B23) 

is obtained for the total system response. 
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APPENDIX C 


ALTERNATE METHOD TO SOLVE THE WIENER-HOPF EQUATION 


References 16 and 17 provide practical methods of solving a given Wiener- 
Hopf equation without finding the explicit form of P(a) when = 0; however, 
those methods are not applicable when x^ > 0. The following is an extension 
to include such a case. 


Preliminary 1 

Let h{s') be a scalar rational polynomial which is strictly proper,^ 


h{s) 


<7(s) 

r(-s) • 1 { b ) 


(Cl) 


where r[-s) = IT (-s + s^.) and 7(s) = IT (s + with > 0, so 

i ^ k ^ 

that s = sr^ and s = -sij^ are poles of h(s) in the RH and LH s-planes, respec- 
tively. Denoting by [ ]^ and [ ]_ those parts of [ ] which are analytic in 
RH and LH s-planes, respectively. Then, when x^ > 0 ^ 


[his^e 


(C2) 


where is a strictly proper rational polynomial in s. The reason is as 

follows. By Fourier *s integral formula, h{s)e~'^9^ is expanded into [ ]^ and 
[ ]_ parts, where 

[h{s)e ^ e dt ^2^^ ^ '^9^^ do 



9vn 




-s + s 




CC3) 


with unique coefficients The Z term just above is denoted as ^(s), which 
has the stated properties. Equation (C2) is valid even when is a multiple 
pole. What is implied by equation (C2) is that, by evaluating [ ] , which is 


lA rational polynomial is proper when the order of the numerator poly- 
nomial is not greater than that of the denominator. If the numerator’s order 
is less than that of denominator, the rational polynomial is strictly proper. 
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a rational polynomial in s, one can avoid difficulty in evaluating the [ ]+ 
part, which is not a simple rational polynomial, but includes terras. 


Preliminary 2 

The matrix ^4(s) of equation (22a) may be modified as follows. 


where 


4(s) = C(-s)yl*(s)C(s) 


/l*(s) 


cf^(-s)cf^(s)C-l(-s)AC-l(s) + f„(-s) • f/(a) 
d^(“s)4z(s) 


(C4) 


(C5) 


Equation (C4) indicates that LH zeroes of det[A{s')'\ are given by 

det[A* = 0, k = I, ku (C6) 

Using equations (C5) and (22b) , one has 

4’^(s)b(s) = C'\s)[A*{s)y^U-s)g{s)Wg^ia) (C7) 

Furthermore, since the matrix i4(s) has its inverse, one has an expression 

4-1 (s) = P-l (s) [P^(-s)]‘l = C'-i(s) [4* (s)]' 1(7-1 (-S) 

(see eqs. (B4) and (C4)). This suggests that p-i(s), and hence v(s) too (see 
eq. (B5)), has a LH factor (7“l(s). 


Possible Expression of v(s) 

In equation (B5) , each entry of vector [P^(-a)]~^b(a) has the properties 
of h{s) in Preliminary 1 above. Hence, putting 

[[P^(-s)]-^b(s)e“V]^ = [p2’(-s)]-ib(s)e‘''^® - ?(s) (C8) 

and substituting into equation (B5) , one has a possible expression for v(s) 
such that 


v(s) = 4 i(s)b(s)e V® - C i(s)p(s) 


(C9) 


where 


(7"l(s)p(s) = P-I(s)c(s) 


(CIO) 



is a rational polynomial vector that is again strictly proper. Following 
Preliminary 2, the LH factor has been artificially inserted in equa- 

tion (C9) . Once the unknown vector p(s) is determined, v(s) is readily solved 
from equation (C9) . Substituting equation (C9) into equation (21) , the W-H 
equation requires 

[z(s)]^ = -[A(s)C"^ (s)p(s)]^ = 0 (Cll) 


Properties of p (s) 

The unknown vector p(s) has the following properties. 

1. In order that v(s) in equation (C9) be analytic in the RH s-plane, 
any singularity of A~^ {s')h{s')e~^9^ in the RH s-plane must vanish with its 
counterpart in 1 (3)p (s) . This implies that RH poles of p(s) must exactly 
coincide with those of i4~^(s)b(s) with the correct multiplicity. 

2. In order that equation (Cll) be valid, that is, z(s) be analytic in 
the LH s-plane, LH poles of C’"^(s)p(s) are not]\ing but LH zeroes of det[A{s')^, 
This fact (ref. 17) is proved as follows. Suppose C“^(s)p(s) has a pole 

s = with Ee{si) > 0. Then after partial fraction expansion, C"^(a)p(a) 
has a term Pi/(a + S]^) with uniquely determined coefficient vector p^^. Exe- 
cuting the [ operation in equation (Cll) gives 

Ad-sO • Pj = 0 (C12) 

which insists, for pi ^ 0, that det[Ad-Si')^ = 0, or equivalently, that 

det[A-^{-Si)] = 0 (C13) 

The statement is then proved by comparing equations (C13) and (C6) . If si 
coincides with a pole of i4(s), a slightly modified proof leads to the same 
result . 


Poles of p(a) 


Properties (1) and (2) determine uniquely the explicit form of p(s); the 
details are presented in the following. When C'(s) and f(s) are given by equa- 
tions (11) and (6), equation (C5) becomes 


A*ds) = 


^ rfin 


frfi 




rrn 


‘^oe '^oe '^a 




(C14) 
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where <^ofie = (s + ]Xf,g), = d^is) , da = dai-s) , etc., are used for short- 

hand notation. Also, 


det[A*{s)'\ 


v^d^i-s) ’ duis) 
dai-s) • dais') 


(CIS) 


where 


* dy^is) - 


^cf ^af ^rfin ■*■ ^Pg^ fnfn (C16) 


and where, with LH zeroes s = -Sy^ (fe = 1, . . ., "ky) , dyis) is so factored 
that 


k. 


u 


duis) = 11 (s + Syj^ , + 2 

k=i 


(C17) 


It must be noted that 


vank[dai-s)da(.s)A*{s)] = 1 ; i = I, 

s - -sa^ 




(CIS) 


which is a consequence of the fact that the output variable y(s) is a scalar 
in this case. 


Again, when C(s) , f(s) and g(s) are given by equations (11) and (6), then 
equation (C7) is given by 


A"l(s)b(s) = C-l(s) 


(■^) dae d^Q /, 


'^of 


9n (®^ 


v^dyC-s)dy(s) 


Wg^is) 


(C19) 


Using equation (19a) and expanding equation (C19) into partial fractions, an 
alternate form 
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CC20) 







ypk \ , Hi 

® ■*■ ^Ukj 3 + dig 


+ 


H2 

(s + 0)^)2 


is obtained with unique coefficients 6^;;, and (either real or complex con- 

jugate), and TJi and 112 (real). Equation (C19) or (C20) indicates that only 
s = k = 1, . . are left as RH poles of 4“^(s)b(s), This is a 

result of hidden pole-zero cancellation. Finally, considering the properties 
of p (s) , one can put 


P(s) 


V ( V - T p- O 

^ V-® * ^Uk s + 


(C21) 


with real (for real or complex conjugate (for complex conjugate 

unknown coefficients and both being two dimensional (/- and 

^-components) vectors. There is a total of 2>^2>^kii unknowns to be determined. 


Unknown Coefficients 3p^ 

Applying property (1) to equation (C9) at s = that is, requiring 

that [v(s) • (-a + ^ with v(a) of equation (C9) and with p(a) of 

equation (C21), one obtains 


= [C^(s)^-l(s)b(s)(-s + (C22) 


When equation (C20) is substituted into equation (C22) , it gives 2xkn unknown 
coefficients 3^^ such as 

^pk = ^ok , k = I, ku (C23) 

which depend upon Xg. 

Unknown Coefficients y^-j^ 

Applying- property (2), equation (Cll) must be evaluated at all LH poles 
of A (s)C“ 1 (s)p (s) , as shown below. At s = k = 1, . . • , ky^ 

[A(s)(7-l(s)]g^_^^^ • Ypk = 0 (C24a) 

at s = -Sa%> ~ 1 > • • •> 
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[(s + Sai')A(.8^C~^(s')pis')] 


(C24b) 


s = 


"®ai 


= 0 


and at 8 


[^(s)C'-l(s)p(s)]^^_ = 0 


(C24c) 


Equations (C24a) , (C24b) , and (C24c) are necessary and sufficient to solve for 
the yp]^. This becomes clear from the following notes. Equation (C24a) is 
simplified to • yp-j^ = 0, but it degenerates into independent 

equations due to equation (C6) . Hence, taking the first row of ^*(s) gives, 
for k = 1, . . . , kj^, 

(^) ~^cf^ofda da^fnfn fn^n\ ' ^pk " ^ (^^Sa) 

^ = ~^uk 

Next, equation (C24b) insists that [(s + s^^)y4*(s)p (s) ] =0, but due to 

equation (C18), it again degenerates into independent relations. Taking 

rp 

again the first row of one has ~ more 

precisely, for i = 1, . . 










" Sai - ^uk 


1=1 


(C25b) 


Finally, equation (C24c) is simplified into 




2 ^oi 

1=1 




(C25c) 


Equations (C25a) , (C25b) , and (C25c) are \l 2 = 2'x-ky^ independent, 

linear and inhomogeneous equations, and determine l^ku unknown coefficients 
Yp^ uniquely with Xg as a parameter. 


Decomposition of v(s) 

Substituting equations (C20) , (C21), and (C23) into equation (C9) , one 
has the solution 


34 



I 


v(s) = (7“^(s) 


IE 


ok ^ "^ok 


^uk ® ^uk 


) ^ 3 - . ^ !!^1 

) ^ * ^g is ^ (0^)2^ 


^ y_pk^\ 

* ^uk ® I 


(C26) 


When equations (C25a) , (C25b) , and (C25c) are solved for = 0, 

~ '^pk^'^g ~ obtained. Using y^i^iO') , equation (C26) is rearranged 

as follows: 


v(s) = g-^Cs) [p^(s)e"’^9'® + Pp(s;x^)] 


(C27) 


where 


p^(s) 


^ ^2 ^ ^ ^ok ~ 

^ * ^g ^ (s + oj^)2 * jf:'i ^ * ^uk 


and 


Pon(^^ 

(S + 0)^)2 


(C28a) 


VpCsiTg) 



e-^9^ - ^ 




s + 


®M?C 



. (C28b) 


Clearly, P^(s) is invariant in x^, and Pn(s;xg') ^ 0 when x~ ->- 0. Also it is 
not difficult to show that p^(s)°is obtained by directly taking the [ ]+ part 
in equation (C8) when x^ = 0, which is the procedure used in references 16 
and 17. ^ 


System Response g^ts) 

Let us consider the system response g-^is) of equation (14) with the feed- 
forward loop closed, which is now decomposed as 


g„(s) = 






(C29) 


where 
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(C30a) 


= 9'C®) - f^(s)p^(s)w^^(s) 

9'wp(®) = -f^(s)Pp(s;T^)w"^(8)e^9'® (C30b) 


Apparently, both have a factor in their denominators; 

these factors are cancelled out with numerators. This is shown as follows. 
Since C'”^(a)p^(s) is a part of v(s) when Tg = 0, and hence satisfies the 
Wiener-Hopf equation (see eq. (21)), 


[^(s)C-l(s)p^(s) - b(s)]^ = 0 


(C31) 


Executing the [ ]+ operation at a = i = 1. 


, reveals that 


[[f^(s)p^(3) - g(s)Wg^is)](_s + 


= 0 


(C32) 




which assures that the factor d^(a) in the common denominator of equa- 
tion (C30a) is cancelled out by its numerator. Thus, one can rearrange equa- 
tion (C30a) as 








CC33) 


where ® polynomial in s, the order of which is not greater than 

that of (s + uig//3) • . Next, gy^is') of equation (C30b) has a factor 


f^(s)Pp(s;T^) 



® ^uk 


■V 


^Pk ^ 

^uk ® ®u/c I 


(C34) 


By virtue of equation (C25b) , it can be shown that equation (C34) is finite at 
s = -S(xi, and i = 1, . . i-Q^, which means that the denominator d^{s) in 

f^(s) is cancelled by the same factor in the numerator of gyp(_s~) . Thus, one 
has a form 
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ku r 

- 

fe=iL 


, „ 1 - eV 

1<} ' ^ok -s + 


^uk 


T “ 


^uk J 


(s) 


(C35) 


Filter Gains 

The following is obtained after straightforward manipulations. As to the 
basic part of the filter, partial fraction expansion of equation (29) with a 
substitution of equation (C28a) results in 



(C36a) 


(C36b) 


^k = 


1 

^9e 



^uk 


0 



1 & — ■ [Vo»: - VpdCO)] (C36c) 

- ^uk 
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where Cq and are real, and is either real or complex conjugate 
same procedure is executed for the additional part of equation (32) , 
one has 












^uk ■*■ 


^uk) 

/3 


2 


• 3 


ok 






X 





^uk 


/3J 


1 - 


= 


^uk 


9 o, 


9e 


^uk 


S - ^uk)^ , , 

^ YpfeCV 


/I 




where Tj^, t-j^ and s-j^ are again either real or complex conjugate. 


The 

and then 


(C37a) 


(C37b) 


(C37c) 


(C37d) 
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(a) Sensor and ride quality measurement stations. 
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(b) Signal flow for extracting gust component. 
Figure 1.- Gust alleviation system. 














^/og, deg/m/sec 


(a) l/\if = 0 (b) = 0.5 sec 

Figure 7.- Effect of lead time on flap control (aft cabin station) 



Figure 8.- Effect of elevator (x^ = 0) . 
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At aft cabin 

— Increases (Xe—0) 
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Figure 9.- Locations of system pole, elevator only. 



w, rad/sec 


Figure 10.- System transfer function 9 'j^Cs) at cabin station 

(t = 0, 1/y =0.2 sec) . 
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(c) Aft cabin station. 
Figure 12.- Concluded. 


E 




a 

o» 



cj, rod/sec 

Figure 13.- System transfer function g A.s) at aft cabin station (elevator only) 
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lncr«o»#$ (X,»0) 



Figure 17.- Locations of system pole when two controllers are used to achieve 

= 0.03 g for Og = 2.1 m/sec. 





Figure 18.- Two 


limiting cases to achieve a given criterion of /n„2 . 
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